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Abstract 

Recently, the interest on theoretical and experimental studies of dynamic properties of the magnetic 
domain wall (MDW) of ferromagnetic thin films with disorder placed in an external magnetic field has 
increased. In order to study global and local measurable observables, we consider the (1 + l)-dimensional 
model introduced by Buceta and Muraca [Physica A 390 (2011) 4192], based on rules of evolution that 
describe the MDW avalanches. From the values of the roughness exponents, global £, local (\ oc , and spectral 
£ s , obtained from the global interface width, hight-diffcrence correlation function and structure function, 
respectively, recent works have concluded that the universality classes should be analyzed in the context of 
the anomalous scaling theory. We show that the model is included in the group of systems with intrinsic 
anomalous scaling (£ ~ 1.5, £ioc = Cs — 0.5), and that the surface of the MDW is multi-affmc. With these 
results, we hope to establish in short term the scaling relations that verify the critical exponents of the model, 
including the dynamic exponent z, the exponents of the distributions of avalanche-size r and -duration a, 
among others. 

1 Introduction 

Due to internal disorder, the dynamic properties of the domains of ferroic materials (ferro-magnetic, -electric, 
or -elastic) is dominated by inhomogeneities. This disorder includes vacancies, defects, impurities and dislo- 
cations, among others. When an external driving field is present, the disorder fixed on the material limits the 
movement of the interfaces or domain walls separating domains. This disorder, quenched in the ferromag- 
netic media, is called Barkhausen noise pQ. Each magnetic domain wall (MDW) exhibits jerky movements 
between multi-metastablc states of quiescence, known as Barkhausen avalanches or jumps. Thermal fluctu- 
ations favor the elastic displacements of the MDW and eventually take it out of a quiescent state, but it can 
only be moved by the external field. Thus, locally the MDW can leave a metastable state by two competitive 
mechanisms: thermal fluctuations or driving fields. 

The study of MDW avalanches has become relevant in the field of memory devices [2H1] > nanowires [SHB] , 
and metal, alloy or semiconductor thin- films [9TH2"]. When ferromagnetic (or ferroelectric) thin- films with 
quenched disorder are placed in an external magnetic (or electric) field, a rich phenomenology is observed 
through the study of dynamic properties in the criticality [T3lfl"5] . An important feature of the avalanches 
in materials with self-organized criticality is their scale invariance in distributions of avalanche-size and - 
duration, with power-law behaviors and well defined critical exponents. These exponents, obtained from 
experiments and models, fall into two classes of universality. One class includes materials where the dipolar 
interactions (or long-range) arc dominant. The other class, in contrast, includes materials where exchange 
interactions (or short-range) are dominant. These properties have been studied intensively in bulk materials 
assuming the usual scaling relationships [T61418J . With the strong interest in the study of ferromagnetic 
thin films, some concepts of 3-dimensional materials started to be reviewed [19] . In thin- films, as well as in 
bulk materials, the movement of the MDW is also dominated by depinning, although the properties start 
to change as a function of the film thickness Experimental results show that below the 200 nm of 
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thickness, the exponents of avalanche-size and -duration have values different from 3-dimensional magnets, 
but also are included in two universality classes based on the dominant type of interactions present in the 
medium [2"lTl2l)] . The theoretical models that adequately reproduce the few available experimental results 
on thin- films, are based on discrete 2-dimensional models [571I2H] ■ However, new issues have arisen when 
establishing relations between different scaling exponents that characterize these systems. 

The pinning-dcpinning transition of these systems is of second order, being the average velocity v of 
the MDW the order parameter. The MDW at non-zero temperature can move taking multi-stable states 
of pinning-dcpinning, i.e. states of quiescence alternated with avalanches, before reaching the phase state 
of moving or pinning. Until today, theoretical approaches to explain the depinning transition of the MDW 
into disordered medium pushed by an external magnetic field, arc based on: (a) the continuous equation 
of Edwards- Wilkinson with quenched noise (QEW) [29H32] and (b) discrete models based on microscopic 
structures and interactions, such as random- field Ising field model with driving (DRFIM) |33fl57] . 

Thin-films and 2-dimensional models can display anomalous behavior in some measured local observ- 
ables, e.g. local interface width or roughness, height-difference correlation functions, and structure function 
or power spectrum, among others. These scaling anomalies in the interface local properties cannot be con- 
cluded directly only from studies of the distributions of avalanche-size and -duration. Recently, a theoretical 
study connected the anomalous roughness exponents to exponents of avalanche, establishing scaling relations 
of general validity, and applying the theory to forced- flow imbibition fronts |38j . Systems with anomalous 
scaling have at least one local observable with nonstationary power- law behaviors for all time [32]. At 
distances or wavelengths much smaller or much larger than the correlation length, this behavior is not nec- 
essarily the same. In contrast, systems with usual or Family- Visek scaling [40] show, in all local observables, 
nonstationary (stationary) power-law behavior for distances or wavelengths much smaller (larger) than the 
correlation length. Each local observable is characterized by at least two exponents. If we study the 2 nd -order 
height-difference correlation or -equivalently- the local interface width, different (equal) exponents of global 
roughness ( and local roughness Cioc lead to anomalous (usual) scaling. Additionally, if we study the struc- 
ture function, also known as power spectrum, different (equal) exponents of global roughness C and spectral 
roughness £ 8 , lead to anomalous (usual) scaling. The roughness exponents £ , Cioc and ( s allow to classify the 
different models and equations that have pinning-dcpinning transitions according to the scaling type. It has 
been shown that the anomaluous roughening can take different forms [4D - PI3] . Ramasco et al. [33] grouped 
the systems into four sets, namely: 



More recent works [33JGH] have confirmed that the DRFIM belongs to a universality class with spatial 
multiscaling and an anomalous scaling that cannot be included in any of the previous categorie^] [45] , since 
C Cioc 7^ Cs and Cioc 7^ 1- On the other hand, the QEW universality class is included in the group that has 
spatial single-scaling and super-rough anomalous scaling [28,46 50]. An interface that evolves according to 
QEW equation is a single- valued elastic string |51!f53] . In contrast, the interface of DRFIM is not single- 
valued as a result of the islands and overhangs left behind by the advancing interface 37,44,54 . Until today, 
we have not identified experimental contributions that report anomalous local properties in the domain wall 
of ferromagnetic thin-films with Barkhausen effect. However, recent experimental contributions have shown 
anomalous scaling in the surface growth of films (e.g. semiconductor and oxides) [55H59| . 

We consider here the evolution model based on rules introduced by Buceta and Muraca [37], and we 
show that it is included in an universality class with intrinsic anomalous scaling. This model takes into 
account the structure and exchange interactions which are present in ferromagnetic thin-films with disorder, 
whose MDW is driven by the external field. Its results for the distribution of avalanche-size and -duration 
are in agreement with experimental data. From the study of the global interface width W we conclude that 
Family- Vicsek scaling is verified at the criticality. We confirm that the pinning-dcpinning transition is clear 
on the saturation of the W. By studying the height-difference correlation function up to 4 th -order, we show 
that the MDW monolayer is multi-affine. As the most important result of this work, we corroborate that the 
model is capable of predicting dynamic anomalies in the local properties of the MDW, that have not been 
experimentally observed but have been predicted theoretically by other models. We show first that the local 
roughness exponent of the height-difference correlation function of second order is different to the global 

First, some authors assumed that the DRFIM belonged to the QEW universality class |28II44| and later, other authors, claimed 
that it had intrinsic anomalous scaling [35] , 
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roughness exponent obtained from the global interface width. This indicates that the Family- Vicsek scaling 
is not appropriate and that an anomalous behavior is present. To complete the analysis, we show that the 
anomalous scaling is intrinsic by determining the spectral roughness exponent, which coincides with the local 
roughness exponent. Finally, we conclude with a discussion about the presented results, our contribution to 
the current state of the subject and the prospects opened by this work. 

2 The model 

We consider a piece of the ferromagnetic thin-film which includes two magnetic domains separated by a MDW. 
On each side of the wall, we assume opposite macroscopic magnetizations in the easy direction. We suppose 
that the medium is composed of magnetic dipoles and randomly distributed point defects, both arranged in 
the nodes of a 2-dimensional square lattice of edge L. In this model we assume that the defects are isolated, 
i.e. two defects cannot be first-neighbors to each other, and the wall in its movement does not include defects. 
The structure of the MDW is considered to be merely a monolayer of dipoles with perpendicular direction 
to the easy direction. The model only takes into account exchange interactions between nearest-neighbor 
dipoles. To simulate the movement of the MDW, the lattice is assumed to have periodic boundary conditions. 
The point defects are represented by a random pinning force rj(i,j), uniformly distributed in [0,1], assigned 
to each node of the lattice. Taking a lattice with density p of dipoles and 1 — p of defects, if r](i,j) < p 
the node has a dipole, otherwise it has a point defect. We characterize the disorder of the lattice using 
the function F(i,j) = 6(p - vihi)), where 0(x) = 1 if x > and @(x) = if x < 0. Since the MDW is 
a monolayer of dipoles, a dipole located at a node (fc, n^) inside the MDW (k = 1, . . . ,L) is described by 
F(k, nk) = 1. A dipole (or defect) outside the MDW is described by F (k, i) = 1 (or 0) with £ ^ rik- 

The evolution rules include only ferromagnetic exchange (or short range) interactions, taking into account 
the balance of the magnetic moment of each side of the MDW in a neighborhood of the point it moves. When 
there is a local unbalance in the opposite direction to the movement, the MDW searches the equilibrium 
with probability c. However, if there is balance, with probability 1 — c, the MDW can only move if there is 
an external force which can remove it from its metastable state. We start the Monte Carlo simulation with 
a flat wall, i.e. initial condition n, = 1 for all i. We introduce the relative position of neighbor nodes with 
respect to the node of the selected column: Xj = n-j+i — rij and yj = rij—i — rij . A randomly chosen node 
(j, rij) of the MDW evolves in the lapse St according to the following rules: 

I. With probability c, if F(j,rij + 1) = 1 and (a) Xj + yj > 2 the selected node is moved one unit, i.e. 
Srij = 1 , or (b) otherwise the selected node is pinned. 

II. With probability 1 — c, if F(j,rij + 1) = 1 and (a) Xj = yj > the selected node is moves one 
unit above its neighbors, i.e. Srij = Xj + 1, or (b) otherwise the selected node moves (or not) to the 
maximum between the neighbors and the same, i.e. Srij = ia&x(xj,yj,0). Also with probability 1 — c, 
if F(j, rij + 1) = and (a) Xj = yj > 2 the selected node is moved one unit above its neighbors, i.e. 
Srij = Xj + 1, or (b) Xj ^ yj and max(xj , j/j) > 3 the chosen node moves to the maximum among its 
neighbors, i.e. Srij = xnax(xj,yj), or (c) otherwise the selected node does not move. 

As the MDW consist of dipoles, a rule is frustrated when a wall point tries to reach nodes with defects. 

3 Scaling analysis and Results 
3.1 Global interface width 

In order to determine the roughness properties of the MDW, we first study the global interface width W as 
a function of lattice size L and time t, defined as 

W(t,L) = {([n j (t)-(n j (t)) L ] 2 ) L } 1/2 , (1) 

where (■ • ■ } = \ S^=i • • • is the spatial average over the system size L and {• • • } is the sample ensemble 
average. Figure [1] (left) shows W as a function of t, for different values of L, close to critical value p c . After 
the early time, the data overlaps much before the crossing time t x = L z with dynamical behavior W ~ t 13 , 
i.e. for t <C t K . In contrast, W ~ const for t 3> t x and the saturation value is function of size L. This growth 
process shows the usual dynamic scaling proposed by Family- Vicsek [40j 

W(t,L) = L<f(t/L z ), (2) 
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Figure 1: (Color online) Left: Plot of the global interface width W(t, L) as a function of time t at several 
lattice size L. The simulation results for c = 0.10 and p = 0.8665 close to critical value are shown. The 
dashed line shows the dynamical behavior W ~ t@ for t <C L z . Right: Plot of the scaling function f(x) 
vs x, obtained via Family- Vicsek scaling [eq. ©] with the left-graph data, taking z = £ = 1.5 . We see 
that the points that do not fit the scaling correspond to the early regime. 



where £ is the global roughness exponent, and z is the dynamical exponent. The scaling function / is defined 
by 

I{X) ~{ const *>1, (3) 

with the growth exponent f3 linked to other exponents by the scaling relation £ = z/3. Figure [T] (right) shows 
the scaling function f(x) and we observe the Family- Vicsek scaling for z = ( = 1.5 . The graphic shows the 
behavior given by eq. (J3]). We notice that the slope of dashed line /3 = 1 is reached with L — > +oo. Figure^] 
shows the change in the saturation behavior of W in the pinning-depinning transition near p c ~ 0.8665. We 
plot W as a function of t, with L fixed, for several values of p close to p c . Below criticality W saturates at 
a constant value, where the average velocity of the wall is equal to zero. In contrast, above criticality W 
has temporal fluctuations around the constant value of saturation, where the average velocity of the wall 
is not zero. Thus for p p c there is a pinning or moving phase, respectively. Figure [3] shows the average 
velocity around the pinning-depinning transition for t t x . The velocity decreases to zero nearby above the 
threshold p* ~ p c as a power law v ~ \p — p*\ e , where 9 is the velocity exponent. Simulations with fixed size 
L and different values of c, let us see that the values of the threshold p and exponent 9 are smooth functions 
of the parameter c. This behavior, characteristic of this model, requires a particularized study that enables 
to compare it with other models. Far above criticality the velocity v ~ p, a characteristic of quenched models 
which should be studied in detail beyond the present work. 

3.2 Height-difference correlation functions 

The movement of MDW, as other surface growth processes, can be studied by the set of height-difference 
correlations functions of m th -order 

G m (r,t) = {(\n j+r (t) - nj (t)\ m ) L } . (4) 

First we consider the correlation G m with fixed time t. We expect a power-law behavior as a function of 
distance r, i.e. G m ~ r m ^ m for r <C £(t) ~ t 1 ^, where £ is the time-dependent correlation length. If Q m 
depends on m the surface is multi-affine. Otherwise, if Cm = C the surface is self-affine. G m is expected to 
saturate for r £. Figure H] shows the correlation G m as a function of distance r at fixed time t for the four 
first values of m. We observe that ( m ~ 1/m for to > 2 and £i ~ 0.84, which ensures that the surface of 
the MDW is multi-affine. The decrease of these roughness exponents with the order to is evidence of spatial 
multiscaling, which usually involves anomalous scaling |59) . In addition, the relation between anomalous 
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Figure 2: (Color online) Plot of the global interface width W(t, L) as a function of time t with lattice 
size L = 2048 for several values of p around the critical value p c ~ 0.8665 and c = 0.10. We notice that 
W for t ^> L z saturates with (without) temporal fluctuations above (below) the critical value. 




Figure 3: (Color online) Log-Log plot of the average velocity v as a function p — p* with lattice size 
L = 2048 and c = 0.10 for p > p* = 0.8666 close to the critical value p c , and temporal data where W 
saturates (t > 10 5 ). We see the characteristic power-law behavior of v close to criticality, i.e. v ~ \p— p* \ 6 
with 9 = 2.01(5). Inside: linear plot of v as a function of p, with the same simulation data. We notice 
that v(p) = for all p < p* . 
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Figure 4: (Color online) Plot of the m th -order height-difference correlation G m (r,t) as a function of 
distance r at time i = 4096, with m < 4. We see the power-law behavior G m ~ r m ^ m for r <C i 1 '*. This 
plot shows results for c = 0.10 and p = 0.8665 close to critical value. The slopes are: Ci = 0.844(4), 
2( 2 = 1.004(4), 3C 3 = 1.036(9), and 4( 4 = 1.045(8). 




scaling and multiscaling has been studied in detail by Asikainen et al. [60l[6T] for single- valued height fronts 
in fractals. On the other hand, there are systems with anomalous scaling and spatial single-scaling [46] ■ Also 
Figure S] shows for r ^> i 1 / 2 saturation of the correlation function G m . Particularly, G2 must have the same 
behavior as W 2 for £ <C r < L (i.e. i « ( x = L z ), since W(L,t) ~ £^ ~ with £ = z/3. We propose a 
scaling 

G 2 M)=^S(r/0, (5) 
where the scaling function Q(x) ~ x 2 ^ 2 if x <C 1 and <?(x) ~ const if x 3> 1. Here £ is the global roughness 
exponent, while C2 := Cioc is known as the local roughness exponent. We should notice the following: 
G 2 ~ r 2 ^ loc £ 2 (C~0°=) for r<(. We have the usual scaling if C = Cioc and, in contrast, we have the anomalous 
scaling if Cioc < (. For r £ there is no difference between the scalings, i.e. G2 ~ i; 2 ^ = t 2< >/ z = t 2/3 . In 
summary, taking £ ~ f 1 / 2 ; we obtain 

G a (r,t)~| t2C/ , r>>e (6) 

Figure O (left) shows G2 as a function of distance r, plotted at different times t. The anomalous behavior 
described here is observed . For r <C £ all curves have given the same power-law behavior as r 2 ^ loc . For 
r £ the curves saturate at different values, as usual. We notice that the crossover points (r = £) of each 
curve has power-law behavior r 2 ^ . Figure [6] shows G2 in the saturation regime as a function of time t. Their 
dynamical behavior is similar to the square interface width W 2 ~ t 2 @ . The growth exponent /3 — 1, obtained 
from the data of Figure [HI is consistent with the value used from the scaling of the global interface width in 
Figure [T] (right). Equation © can be obtained from eq. ((5} or using the scaling 

G 2 (r,t) = r 2 <g A (rA 1 /*), (7) 

where the so-called anomalous scaling function is 

1 \ / U ~ K1 u<^l , . 

gA ^ ~ 1 Ol, (8) 

being 

Kl = 2 (C - Cioc) (9) 

k 2 = 2C- (10) 

Figure [5] (right) shows the anomalous scaling function gA = r _2< =G'2 as a function of u = rt~ x l z . Taking 
z = C = 1-5 , the scaling is very good and confirms the anomalous feature of the model. We also find that 
Cioc 0.5. 
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Figure 5: (Color online) Left: 2 nd -order height-difference correlation function G%{r, t) vs. distance r 
at different times t. The simulation results for c = 0.10 and p = 0.8665 close to critical value can be 
observed. We take the lattice edge L = 1024. For r <C £, we observe that G2 ~ r 2< » loc and that the curves 
are parallel between them, which is characteristic of anomalous behavior. From eq. ([7J for r = £ we 
conclude that G2 ~ r 2l » as this plot shows. Right: Plot of the anomalous scaling function gA{u) = r 2 ^ G2 
as a function of u = t l / z r, obtained from scaling of the data used in left-graph at various times. The 
exponents use for scaling arc ft = 1 and £ = 1.5. The obtained slopes (— «i and — K2) are in agreement 
with this election. With these values we obtain (\ oc = 0.495(20). 
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Figure 6: (Color online) The plot shows the 2 nd -order height-difference correlation G2 as a function of 
time t for £ <C r < L [saturation regime of Figure 3 (left)]. Results for c = 0.10 and values of p close to 
estimate p c = 0.8665(5) can be observed. The data follows a power-law similar to the one which follows 
W 2 for t <C L z . The power-law fit allows to determine (3 = 0.996(39). 
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Figure 7: (Color online) Left: Plot of the spectral function S(k,t) vs. wave number k at various times 
t. There is a clear power-law behavior and an anomalous character because the curves do not overlap 
in this regime. Right: Plot of the spectral scaling function ^(u) — k 2 ^ +1 S as a function of u = t x ^ z k, 
obtained from scaling the data used in left-graph. To perform the scaling, ( = 1.517 and /3 = 0.992 
were taken. The obtained slopes {k% + 1 and K3) are in agreement with the selected values of £ and /3. 
Accordingly, ( s = 0.51(3) has been obtain. 



3.3 Structure function 

Instead of considering the autocorrelation function of the height, it is convenient to introduce its Fourier 
transform, known as the structure function or spectral power density S(k,t). Introducing the Fourier trans- 
form of wall points hk(t) — J2jLi n j(t) e lfcj 'j with wave number k = 2ni /L (i = 0, 1, . . . , L — 1), the surface 
structure function is 

S(M) = (n fc (i)n_ fc (t)) (11) 

Figure [7] (left) shows the structure function S as a function of the wave number k at different times t. We 
observe that for intermediate wave numbers there is a clear power-law behavior and non-overlapping curves. 
This is characteristic of anomalous scaling functions. It is appropriate to propose that the structure function 
scale as 

S(k,t) = k-( 2i+d) V(t 1 / z k) (12) 
where the so-called spectral scaling function is 

^-\ U T U<< ] , (13) 

with d the surface dimension |43| and 

k 3 = 2(C - C.) • (14) 

Figure [7] (right) shows the anomalous spectral scaling function ^ = fc 2 ^ +1 5 as a function of u = kt x / z . The 
anomalous scaling fits very well and £ s ~ 0.5 has been determined, which justifies its intrinsic character. 

In all figures, c = 0.10 has been chosen arbitrarily, although other values of c G (0, 1) show qualitatively 
the same results. Furthermore, the critical value p c is a very smooth function of c. 



4 Conclusions 

The main contribution of this study was to characterize the statistical properties of local observables of 
MDW, such as the correlation and the structure function. We use a discrete rule-based model to describe 
the motion of the MDW between two domains of a ferromagnetic thin-film sample with defects, when placed 
in an external magnetic field. This 2-dimensional model has successfully described MDW avalanches [2"7] 
and has shown a pinning-depinnnig transition whose properties had not been characterized before. With the 
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purpose of obtaining global properties we considered, as usual, the time evolution of the interface width for 
different lattice sizes. We observed Family- Vicsek scaling at criticality with z = ( = 1.5 . Around criticality, 
the global interface width in the saturation displayed the characteristic of the pinning-depinnig transition. 
Above (below) criticality the global interface width saturated with (without) fluctuations. We found that 
the interface is multi-afSne by studying the height-difference correlation functions of different orders, all 
with power-law behaviors. We show that the 2 nd -order correlation exponent (or local roughness exponent) 
is not equal to the global roughness exponent. This fact, among others outlined here, made us conclude that 
the scaling is anomalous. We performed the scaling with z — £ = 1.5 to determine (i oc — 0.5. With the 
purpose of classifying the anomalous scaling, we studied the structure function. In the regime of intermediate 
wave number, the power-law scaling was clear and allowed us to determine, with high precision, the spectral 
exponent and ensured that £ s — (j oc . Following the criteria introduced by Ramasco et al. [43], we were able to 
affirm that our model belongs to the set of systems that have intrinsic anomalous scaling. We have carried out 
a rigorous analysis to ensure the existence of intrinsic anomalous scaling. The local and spectral roughness 
exponents, obtained from regressions, coincide within error range determined by the errors of the calculated 
exponents. If we analyzed the rules of evolution, in detail, we would observe that the MDW is formed by 
terraces and plateaus. This is a feature of interfaces with intrinsic anomalous roughening. Systems with 
interfaces that leave islands and overhangs in their advance as the DRFIM are not included in this category. 
In the context of the study of the Barkhausen effect in ferromagnetic thin-films, we believe it is promising 
to establish the connection between roughness and avalanche exponents, following the known theoretical 
framework |38j . According to dynamic rcnormalization group results, the intrinsic anomalous scaling cannot 
occur in homogeneous (i.e. non-disordered) and local growth models |62j . In agreement with this result, the 
quenched disorder in our system is responsible for intrinsic anomalous roughening. The situation is similar 
to what occurs in other systems with a depinning transition, as fluid imbibition in disordered media. Finally, 
another outstanding point is the morphology of the domain wall from the set of local and global exponents, 
which should be subject of future studies for models and experiments. 
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